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Abstract 

Glass-ceramics with nanocrystals present a transparency higher than that expected from the 
theory of Rayleigh scattering. This ultra-transparency is attributed to the spatial correlation 
of the nanoparticles. The structure factor is calculated for a simple model system, the random 
sequential addition of equal spheres, at different volume filling factor. The spatial correlation given 
by the constraint that particles cannot superimpose produces a diffraction peak with a low S(q) 
in its low-q tail, which is relevant for light scattering. The physical mechanism producing high 
transparency in glass-ceramics is demonstrated to be the low density fluctuation in the number of 
scatterers. 
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Glass-ceramics with nanometric size of the crystals are of great interest for photonics, 
since they combine the spectroscopic properties of crystals with the mechanical properties 
of glasses [l, 2, [j, [l]]. Many glass-ceramics, in particular those obtained from oxyfluoride 
glasses, have shown a transparency similar to that of glasses, much higher than that expected 
from the theory of Rayleigh scattering . Recently, Edgar et al. measured the optical 
extinction coefficient for fluorozirconate glasses containing a volume fraction of about 0.2 
of barium chloride nanoparticles, finding that the glass ceramics are about a factor of six 
more transparent than predicted by the Rayleigh scattering theory, based on estimates 
of particle sizes from X-ray diffraction and electron microscopy measurements jfj. They 
modelled the electric response of the material within a discrete dipole approximation, which 
accounts for the particle-particle interaction. However, they did not succeed in reproducing 
the observed quenching of the attenuation. In our opinion the dipole-induced-dipole effects 
are not crucial for the scattering intensity. The particles, in a first approximation, can 
be considered independent and excited by a homogeneous external field. The main effect is 
related to their spatial arrangement and the consequent interference pattern of the scattered 
field. Tick suggested that the low attenuation is related to this effect, but quantitative 
calculation are lacking. Hopper treated the case of a spinodal decomposition and Hendy 
found a g 8 ~ A -8 dependence of the Rayleigh cross section in the low q limit, q being the 
wavevector of the light (q = 2tt/\) 8|. 

In the absence of any spatial correlation the scattering attenuation is given by 



a 



where the sum in on all the particles in the volume unit and Oi is the cross section of 
Mie scattering of the i-th particle. If the interaction among particles is neglected, i.e. each 
particle is considered as surrounded by a homogeneous medium, and if all particles are equal, 
Eq (1) reduces to a = pn&, where pn is the particle density. Furthermore, if the size of the 
particle is much smaller than the wavelength of the light, Rayleigh scattering occurs. The 
cross section for Rayleigh scattering by a not absorbing spherical particle of radius R and 
refractive index n n , embedded in a not absorbing medium of refractive index n m , is given 
by H: 
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where m = n p /n m , k = 2ix/X = 2irn m /Xo is the wave vector of the light of wavelength A in 
the medium and Ao in vacuum. 

In the presence of spatial correlation, interference effects occur and the total scattered 
field at any angle (9) of scattering, or at any exchanged wavevector q = (An/X) sin(#/2), is 
the sum of all individual fields. For N equal particles excited by a plane wave, the individual 
scattered fields are given by: 

Ei(q) = E (q)exp(zq ri ), (3) 
and the scattered intensity is given by: 

E(q) 2 = iVE (q) 2 %), (4) 
where S(q) is the structure factor of the system: 

%) =| ^exp^qr,) \ 2 /N. (5) 

i 

We will show that any spatial correlation lowers significantly S(q) at the small q- values 
with respect to the value S(q) = 1 of the uncorrelated system. The attenuation of the 
light will depend on the scattering at all angles (0 < 9 < n) or q- values (0 < q < An/X). 
Furthermore, we will show that for small q- values, i.e. the ones which are relevant to Rayleigh 
scattering, the structure factor is nearly constant, i.e. S(q) ~ S(q — 0) = So- Therefore the 
attenuation, in the presence of a particle correlation measured by S(q) is well approximated 
by a c = a s S = PN^RaySo- 

An unavoidable correlation among equal spherical hard particles is given by the impos- 
sibility of interpenetration. Interestingly, this effect may enforce non-trivial spatial corre- 
lations even on length-scales larger than 2R. A simple process yielding three dimensional 
systems subjected only to such hard-sphere constraint is given by the random sequential 

r~i 

addition (RSA) [lOj]. In such process, the n-th particle is placed randomly within a cubic 
box with periodic boundary condition only if it finds the space left free from the previ- 
ous n — 1 particles. The correlation among the particles increases when the filling factor 
cf) = A/3ttR 3 pn, where R is the particle radius, grows. For = 0, there is no spatial 
correlation, while the maximum filling factor for RSA in 3-D is 4> sat — 0.3828 [lo1 |. 

Figure 1 shows the calculated S(q) for various values of 0, obtained by averaging over 
30 realizations of the RSA process with N = 10000 particles. As <fi increases, a diffraction 
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peak, located at q p ~ 2ir/2R becomes more and more well defined. The radial distribution 
function g(r) (not shown here) is zero for r < 2R and presents a single peak at r p > 2R. 
The radial correlation decays very fast with r for r > 2R for any 0, similarly to the case 
near (f) sat (0 = 0.3812), described by Torquato et al. [lo| . 
Figure 1 around here 

The quenching of the S(q) in the low-g region is shown in more detail in the inset of Fig. 
1. When plotted as a function of q 2 , the parabolic behavior [h]] 

S(q) = S + aq 2 (6) 

appears as a straight line. Here S is the low-g limit, being q min = 2n/L, where L is the edge 
of the cubic box, the minimum q- value accessible in the calculations. The fitted So values 



are reported in Fig. 2, together with the value found near the saturation density 10]. 

A quantitative prediction can be obtained by studying the simplest form for g(r) which 
takes into account both the hard-sphere constraint and some possible correlation over a 
lenght-scale £, 

g (r) = r <2R 

g{r) = 1 + Ae~ r/ t r > 2R. (7) 
In fact with that choice one finds: 

So = 1 - 80, (8) 

which suggests that the increase of the filling factor is the main factor which determines 
the drop of Sq. Therefore, the very short range correlation of our model system, which 
only avoids particle superposition, is able to reduce the scattering intensity by about an 
order of magnitude for in the range of 0.2 — 0.3, values that are typical of glass-ceramics. 
Furthermore, no medium or long range correlation seems to be needed for quenching Rayleigh 
scattering. But, in fact, some degree of long range order is produced. A deeper insight is 
obtained by the following study. The box containing N particles (N= 10000 in our case) can 
be decomposed in n cubic smaller boxes, each containing M = N/n particles, on average. In 
the absence of any particle correlation, a Poisson distribution in the number Mj of particles 
will be present. For large values of M, it is well approximated by a Gaussian distribution 
with variance o~ 2 M = M. In the inset of Fig. 2 we show the results for n = 125, M = 80 
and for two values of the volume filling factors, <fi = 0.05 and <ft — 0.3. The variance of the 



particle number is smaller than M and decreases as the filling volume increases. This can be 
easily understood, since in the RSA method of box filling the n-th particle will have larger 
probability of being accommodated in a box with smaller particle density. This will produce 
an equalization of density, reducing its fluctuations over large scales to lower values than 
those that would be present for a random distribution. Figure 2 summarizes the results of 
the calculations for different n-values of the box partition, n = 125, 64, 27 and 8, showing 
the normalized variance a' M 2 = <7m 2 /M. 
Figure 2 around here 

The above feature is not limited to RSA systems. As a matter of fact, in a system of 
classical particles, S is always related to the fluctuations of the number of particles M 



within spheres of (large) radius / by the relation 111 ]: 



(M 2 ) - (Mf 



1 + Aixp N / drr 2 (g(r) - 1) 
Jo 



(M) 

- S , (9) 

where the limiting value is approached when I — > oo (this explains the observed size depen- 
dence of <j' m on the sampling size M in Fig. 2). 

In real systems, as also suggested by the model of Eq. 7, So can be much lower than what 



is found in RSA 



Is 



A naive argument, which basically exploits Eq. 9, allows accounting 
for such a decrease. In fact, even though the processes of nucleation and coarsening are 
very system dependent, the growth of nanosized particles usually occurs by diffusion within 
limited spatial region, each particle having its own basin for gathering material. In a volume 
with a characteristic distance of the order of 1/im relevant for Rayleigh scattering, we will 
find a certain number N BB of the building blocks (BB) that will compose the nanoparticles. 
If we assume that M = N BB /n BB identical nanoparticles, each one containing n BB BB, 
are formed by limited diffusion within the given volume, M will have the same relative 
fluctuation as N BB : <Jm/M = &n bb /N bb . The fluctuations on the particle number will 
therefore be proportional to l/n BB and, by Eq. 9, so will be So and the intensity of the 
scattered light. Note that since the cross section a Bay is proportional to the square of the 
particle volume (see the R 6 dependence in Eq. 2), o~R a ypN is proportional to the particle 
volume (assuming that all the BB have been precipitated in particles and n BB M = Nbb), 
and finally a c = CRayP^So becomes size independent. Therefore, no increase of attenuation 
is expected with the increase of the particle size, since the increase of the cross section for 



Rayleigh scattering is exactly compensated by the decrease of the structure factor. 

This simple model does not account for other important effects, as the size distribution 
of particles or the formation of more ordered structures. In any case, the physical mecha- 
nism producing high transparency in glass-ceramics even for relatively high particle size is 
demonstrated to be the low fluctuation in the number of scatterers, due to spatially limited 
diffusion during coarsening. 

The above results indicate two connected methods for estimating the quenching of the 
Rayleigh scattering in a glass-ceramics. The structure factor can be measured by Small 
Angle X-ray scattering. If the measurements are extended, as far as possible, down to low-q 
side of the peak, the very \ow-q S(q), So, should be possibly obtained by extrapolation to 
the parabolic behavior of Eq. 6. The fluctuations in the spatial particle distribution can be 
measured by counting the particles using TEM tomography. The decrease of Mie scattering is 
expected to be reduced by the factor So = a' M 2 with respect to that calculated for a random 
distribution of particles. The measure of light attenuation and of the structure factor, 
combined with the calculation of the Mie cross section for the single particle scattering, 
should quantitatively confirm this explanation of the observed ultra-transparency of glass- 
ceramics with nanocrystals. 

The authors acknowledge M. Ferrari for useful discussions about glass-ceramic materials. 
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FIG. 1: Calculated structure factor for 3-D random sequential addition of equal spherical particles 
and for different volume filling ratios, (j). The arrows indicate the values of q = 2tt/2R, where R is 
the radius of the spheres. The inset shows the low-q part as a function of q 2 . 
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FIG. 2: Calculated low-q limit of S(q) (dots). The star shows the value (Sq =.05) at the saturation 
density from ref. [h]]. The dashed line shows the ^-dependence given by Eq. 8. Normalized 
variance of the particle density in cubic boxes containing 10000/8 (open squares), 10000/27 (circles), 
10000/64 (triangles) and 10000/125 (inverted triangles), on average. Inset: distribution of number 
of particles in the 125 small cubes of size L/5 for volume filling ratios 0=0.05 and 0.3, calculated 
over 48 and 27 samples, respectively. The curves are fits by Gaussian curves. 
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